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ABSTRACT 

A computer simulation of a flywheel energy 
storage single axis attitude control system is 
described. The simulation models hardware which 
will be experimentally tested in the future. This 
hardware consists of two counter rotating flywheels 
mounted to an airtable. The airtable allows one axis 
of rotational motion. An inertia DC bus coordinator is 
set forth that allows the two control problems, bus 
regulation and attitude control, to be separated. 
Simulation results are presented with a previously 
derived flywheel bus regulator (Kascak, 2001) and a 
simple PID attitude controller. 


INTRODUCTION 

Flywheel system and technology development 
work has been an ongoing activity at NASA Glenn 



FIGURE 1 - HIGH ENERGY FLYWHEEL 
FACILITY 


Research Center over the past several years. Most 
recently, work has been conducted under the 
Aerospace Flywheel Technology Program (AFT) and 
the Flywheel Energy Storage System Program 
(FESS). This paper discusses simulation results of a 
control algorithm for combined power bus regulation 
and single degree of freedom torque and momentum 
control using a two flywheel system. A discussion of 
related programs sponsored by NASA and the US Air 
Force has been done. (McLallin, 2001). 

The algorithm presented will be experimentally 
verified in the Attitude Control and Energy Storage 
Experiment Phase 0 (ACESE_0). In this experiment 
two technology demonstrator flywheel modules 
developed under FESS and a previous flight program 
will be integrated at NASA GRC’s High Energy 
Flywheel Facility (HEFF). Figure 1 shows the facility 
which consists of an airtable where the modules are 



Air Table 

FIGURE 2 - HIGH ENERGY FLYWHEEL 
FACILITY SCHEMATIC 
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mounted surrounded by a water containment safety 
system (Trase 2002). Thermal, vacuum, and 
instrumentation support hardware can be seen on the 
upper platform. A simple schematic of the setup is 
shown in Figure 2. The experiment will be a 
hardware demonstration of a flywheel system which 
provides both power bus regulation and a single axis 
torque and attitude control. The long term objective is 
to extend this work to a bus regulation and three 
degree of freedom attitude control system 
representative of a satellite platform. 

NOMENCLATURE 

I DC - DC current provided by flywheel system (A) 

I DC m - DC current into flywheel inverter m (A) 

I p/s - DC current out of the power supply (A) 

J - moment of inertia about axis x (kg*m 2 ) 

L d - direct axis self inductance (H) 

L q - quadrature axis self inductance (H) 
p - derivative operator 
P - number of poles 
T d - disturbance torque (N*m) 

T - electromagnetic torque (N*m) 

T p - platform torque (N*m) 

V DC - DC bus Voltage (V) 

- mutual flux linkages by rotor magnets, (V*s) 
C0 ; 2 - electrical frequency (rad/s) 


SYSTEM MODEL 

The ACESE_0 simulation is divided into two 
parts, the controllers and the plant models (Figure 3). 
Each controller regulates progressively higher level 
states as we move from the inner to outer control 
loop. Plant models are provided for each hardware 
subcomponent. The model does not include 
switching effects, motor torque ripple, or other higher 
order effects that are required to estimate power or 
torque quality. The purpose of the model is to 
develop a functional controller for the hardware 
demonstration. 

Plant models were developed for the two flywheel 
modules, a single degree of freedom platform 
mounted on an airtable, a power supply, an electrical 
load, and the DC bus. 

The controller has two nested loops and a 
coordinator. The innermost loop is a motor current 
regulator for each flywheel. The outer loop has a DC 
bus regulation algorithm and attitude control algorithm 
operating simultaneously. The system level 
commands of bus voltage and pointing angle get 
issued to the highest level loop. The coordinator 
parses the demands from the bus regulator and 
attitude controller into current commands for each 
flywheel module. 

Inputs to the system controllers include: V* DC , 
f charge (to DC Bus regulator) and theta* (to attitude 
regulator). The parameters that are varied are: V com , 
I/imit (power supply plant model), load resistance (load 
plant model), and disturbance torque (platform plant 
model). 

Mechanical Plant Models 

The flywheel modules which will be incorporated 
into HEFF are called D1 and HSS. The mechanical 
plant models included in the system are the D1 
flywheel module, the HSS flywheel module, and the 
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single degree of freedom plant model. 

The D1 Flywheel Module (D1) and the High 
Speed Shaft (HSS) flywheel module used in this test 
are not identical, however they both have generally 
the same topology. They each have magnetically 
suspended rotors with permanent magnet 
motor/generators, mechanical touchdown bearings, 
and an integral housing/vacuum enclosure. Table 1 
contains the flywheel module properties. 

The mechanical system is modeled as a single 
degree of freedom second order system. The airtable 
and all the components mounted on it other than the 
flywheel rotors will be referred to as the platform. The 
platform components are modeled as inertias about 
the spin axis of air table. These components are the 
base, avionics, and flywheel stator assemblies. The 
flywheel rotors are included as torques acting on the 
platform. Figure 4 shows a free body diagram of the 
system. Axes 0,1, and 2 represent the rotation axes 
of the platform, flywheel 1, and flywheel 2 
respectively. The equation of motion for the platform 
is: 

J0 O +C0 O +kQ 0 = (1) 

The work in this paper does not include platform 
damping and stiffness effects; these are set to zero in 
the model. The terms are included for future 
modeling of the cable harness stiffness and bearing 
damping in the experimental hardware. 

The total inertia about the 0 axis, J 0 , is the sum of 
the inertia of the base, avionics, flywheel stator 
assemblies and flywheel rotors. 

Jo = J a +J b + 

(J, 1 + m r\ X l ) + ( J , 1 + m si X l ) (2) 

(J r2 m r 2 X 2 ) (J s2 m s2 X 2 ) 

The base and avionics are modeled as circular 
cylinders. The flywheel stators are considered 
circular cylindrical shells located on a parallel axis. 
The flywheel rotors are circular cylinders on parallel 
axis. The standard equations are applied (Meriam 
and Kraige, 1986). 

J zz = j tnr 2 circular cylinder (3) 



FIGURE 4 - FREE BODY DIAGRAM 
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J zz =mr 2 circular cylindrical shell (4) 

J axis_a = J axis_b + ™ d a, Jo _b P arallel aX ' S ( 5 ) 

The torques applied to the platform are equal and 
opposite to the torque applied to the rotors of each 
flywheel model. Resultant torques from rotor 
translation are neglected. The platform inertia reacts 
to the sum of the torques. 

Y J T=T p =-T eX -T e2 (6) 

The flywheel module is modeled as a single 
degree of freedom rotor with a torque applied by the 
motor/generator. The motor is modeled as a 
synchronous permanent magnet machine (Krishnan, 
1996). 

Electrical Plant Models 

The electrical model includes: the DC bus voltage 
plant model, the power supply plant model, and the 
load plant model. 

The DC bus is modeled as a capacitor with 
various current sources charging or discharging it 

TABLE 1 - SYSTEM PROPERTIES 


Value Property 

System 

Jo 

25.22 kgm 2 

platform inertia about axis 0 

c 

0 

platform damping 

k 

0 

platform stiffness 

D1 Flywheel Module | 

m r1 

17.2 kq 

rotor mass 

Jn 

0.0664 kg m 2 

rotor polar moment of inertia 

m s i 

61.2 kg 

stator mass 

Jsl 

560 kg m 2 

stator polar moment of 
inertia 

Xl 

0.254 m 

distance from axis 0 to 1 

Lql 

0.000025 H 

q axis inductance of motor 

L d , 

0 .000019 H 

d axis inductance of motor 

Rsi 

0.0145 Q 

l-n stator resistance of motor 

^-afl 

0.00733 V-s 

rotor magnet flux linkage 

Bi 

0 

viscous drag on rotor 

| HSS Flywheel Module | 

m r2 

6.53kg 

rotor mass 

Jr2 

0.00377 kg m 2 

rotor polar moment of inertia 

m s2 

56.7 kg 

stator mass 

Js2 

518 kg m 2 

stator polar moment of 
inertia 

x 2 

0.254 m 

distance from axis 0 to 1 

L q2 

0.000142 H 

q axis inductance of motor 

L q2 

0.000101 H 

d axis inductance of motor 

Rs2 

0.035 Q 

l-n stator resistance of motor 

^-af2 

0.014456 V-s 

rotor magnet flux linkage 

b 2 

0 

viscous drag on rotor 

Base | 

m b 

40.1kg 

base mass 

Jb 

2939 kg m 2 

base polar moment of inertia 

Avionics | 

m a 

40.1kg 

avionics mass 

Ja 

2939 kg m 2 

avionics polar moment of 
inertia 

C 

5860 pF 

dc bus capacitance 






(Figure 5). These currents are the inputs to the DC 
bus model and consist of: the DC currents provided 
by the flywheels ( I DCI + the power supply current 
(Ip/s ), and the load current (I , oad ). 



FIGURE 5 - DC BUS MODEL 

The power supply is modeled as a voltage 
regulating current source. The power supply model 
has two inputs, one for a voltage setpoint and one for 
a current limit. The power supply will provide 
whatever current, I P/S , is necessary to keep the bus 
voltage at the voltage setpoint as long as the current 
is less than the current limit. To accomplish this in the 
model the error between the voltage setpoint and the 
bus voltage is taken. This error is run through a PI 
control law to provide a current command. This 
current command is compared to the current limit. If 
the command is less than the current limit, I P/S is equal 
to the current command. If the current command is 
greater than the current limit, I P/S is equal to the 
current limit. 

The load is modeled as a resistor. Thus the 
current consumed by the load is just the DC bus 
voltage divided by the value of the resistor. 

Controller Model 

The controller model includes: a current regulator 
for each flywheel, a DC bus regulator, an attitude 
regulator, and an inertia and DC bus coordinator 
(Figure3). 

The two motors are controlled by field orientation, 
with back emf voltage decoupling (Krishnan 1999). 
The d-axis rotational reference frame current is 
commanded to zero and the q-axis rotational 
reference frame current command varies to provide 
attitude control and bus regulation. A detailed 
description of the innerloop controller, albeit without 
voltage decoupling, can be found in Kenny (2001). 

There are three modes of DC bus regulation: 
charge, charge reduction, and discharge. The charge 
mode represents operation in full sun when the 
photovoltaics system provides current and voltage 
regulation. In the charge mode the DC bus regulator 
charges the flywheel system at the charge setpoint 
(/ ch arge)- The charge reduction mode represents 

operation in partial sun when the current generated by 
the photovoltaics system is not sufficient to meet both 
the load demand and the flywheel charge setpoint. 
During this mode of operation the flywheel system 
regulates the DC bus voltage to V ic . The discharge 
mode represents operation in eclipse when the 
photovoltaics system provides less current than 
required by the load. In both charge reduction mode 
and discharge mode the DC bus regulation algorithm 
uses the DC current from each of the two flywheel 
inverters to regulate the bus voltage. More detail on 


the DC bus regulator and the modes of operation can 
be found in Kascak (2001). 

The q-axis rotational reference frame current 
commands to the two flywheel motor controllers, C\ sl 
and i* f q s2 > provide both attitude control and DC bus 
regulation. The first relationship necessary to 
determine the appropriate current commands is the 
power balance. The AC power in one of the flywheel 
motors can be expressed as shown in Novotny et. al 
(1996). 


^ ^ / i* • r r * i* \ 

flywheel ^ ' qs^ qs ^ds^ds ' 


Exploiting the fact that the d-axis current is 
regulated to zero and equating the AC and DC power, 
the following is obtained where I DCm is the DC current 
into one of the flywheel motor inverters. 

= |(v;c> (») 


The q-axis voltage equation is (Krishnan, 1999) 


v ' qs =iqs R s + L qsP i qs + *A°>rA& + “ A/ ^ 

In steady state the derivative term goes to zero. Also 
the d-axis current is regulated to zero. Finally the 
resistive term is much smaller than the back emf term 
especially at high speeds, thus we have (Kascak 
2001), 

R DC^DC,m ~ ^ (® r^afiqs ) ( 


Equation (11) expresses the DC current in terms 
of the rotational reference frame q-axis current, f qs , 
with a dependence on a motor constant, X af the rotor 
speed, and the DC bus voltage V DC . 


3co ,X 


af 


1 DC,m 


2V, 


( 11 ) 


DC 


The total DC bus current, I DC , is the sum of the 
two motor inverter DC currents, I DCI and I DC2 . This 
can be expressed as shown in (12) in terms of the 
rotational reference frame q-axis current, f qs , for each 
motor. 


I DC ~ c \s\ + ( 12 ) 


qs 2 


3o 


c = 


A/1 


2V, 


d = 


3co, 


DC 

A a/2 


(13) 


2V n 


DC ( 14 ) 

The torque applied to the platform is equal and 
opposite to the torque applied to the rotor. The 
expression for torque produced by the flywheel PM 
motors is given as (Krishnan, 1999) : 

3P\ 


e 2 2 L 


af 


+ 


(L d -L q y ds }•; (is) 
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However the rotational reference frame d-axis current 
is controlled to be zero. Thus the torque can be 
expressed as: 

T e = ffV'i ( 16 > 

The total torque that is applied to the platform is the 
negative sum of the two motor torques, T el and T e2 - 
This is expressed (from 6) as: 

T p = ~ T el ~ T e2 = -«Cl “ bi 'qs2 ( 1 7 ) 


2 2 
3^ 
2 2 




af\ 


X 


af 2 


(18) 

(19) 


By solving equations (12) and (17) 
simultaneously, the control variables i r qsI and i r qs2 can 
be found which provide the desired DC bus current, 
I DC , and the platform torque, T p . The result is given in 
equations (20) and (21). 


l 


qs\ 


bI dc + dT p 
cb-da 


^qs2 


-b 


( 20 ) 

( 21 ) 


Thus the required values of i' qsI and i' qs2 to provide 
the desired amount of DC current and the desired 
amount of torque to the platform are found. 


RESULTS 

System model simulations were run for three 
different test cases: DC bus regulation throughout 
mode transitions with a constant platform angle, 
platform angle and DC bus response to a step change 
in commanded angle with a constant electrical load, 
and the system response to a step change in 
disturbance torque. Commanded variable and 
system parameter settings for each of these 
simulations are summarized in Table 2. 


DC Bus Regulation 

This simulation verifies that the controller maintains 
DC bus voltage regulation and platform angle 
regulation as the system transitions through all three 
modes of operation: charge, charge reduction, and 


TABLE 2 - SIMULATION PARAMETERS 



Commanded Values 

Bus Condition 

Test 

Bus 

Regulator 

Attitude 

Regulator 

Power 

Supply 

Load 

DC Bus 
Regulation 

V' DC = 130 

I* charge^ 5.4 

theta’=0 
Tj= 0 

6-0 kW 
ramp 

2 kW 

Angle 

Regulation 

V' DC = 130 

I*charge=^A 

theta’= 

1 .0° step 
T d =0 

0 kW 

0 

kW 

Torque 

Disturbance 

Rejection 

V' DC = 130 

/W' 5 - 4 

theta’=0 
T d =0.05 
N-m step 

0 kW 

0 

kW 


discharge. The mode transitions were achieved by 
starting the model in charge mode with the power 
supply current limit set to maximum (from t=0 to t=1 
sec.), then decreasing the limit linearly to zero (from 
t=1 to t=3 sec.). This represents the spacecraft 
moving out of full sun and into eclipse. 

While examining bus performance, it is helpful to 
refer to (22), which is derived from Figure 5: 

I load = Ip/S - IdCI — IdC2 (22) 

Power modes for this simulation are displayed in 
Figure 6. For the first second of the simulation, the 
system is in charge mode; the power supply regulates 
the bus voltage. The do bus voltage is held at 135 
volts by the power supply. The two flywheel units are 
being charged at I DCI and I DC2 respectively where I DCI 

^ IdC 2 I charge- 

At t=1 second, the system transitions into charge 
reduction mode as the power supply current limit is 
reduced. At this point the sum of the load current, 
I had , and the commanded charge current, f charge , is 
greater than the amount of current the power supply 
can provide and the power supply can no longer 
regulate the bus voltage. The flywheel system takes 
over bus regulation, as can be seen by the drop in V DC 


DC Bus Regulation Simulation: DC Bus Response 



time (sec) 


FIGURE 6- DC BUS RESPONSE DURING 
REGULATION 



time (sec) 


FIGURE 7- PLATFORM ANGLE DURING DC 
BUS REGULATION 
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to the flywheel system voltage regulation set point 
value. This setpoint, V DC , is 130 volts but due to the 
large capacitor in the model it takes approximately 0.1 
seconds for the system to reach its final value. 

The system transitions from charge reduction to 
discharge mode when the power supply current can 
no longer meet the load current requirements and I Da 
and I DC2 swing negative. At this point, the flywheels 
are discharging; current from the flywheels is required 
to meet the needs of the load. In this mode, the 
flywheel system is still maintaining DC bus regulation. 
In this simulation, this transition occurs at about two 
seconds. In the modes when the flywheel regulates 
the DC bus the voltage is held within 0.1 volts of the 
setpoint. 

Figure 7 displays the platform angle change 
during the mode transitions. There is no change in 
the commanded or actual platform angle during the 
mode transitions. 

Angle Regulation 

The angle regulation simulation examines 
platform angle and bus voltage regulation as the 
system responds to a step change in commanded 
platform angle. This represents a commanded 
spacecraft direction change. For this simulation, the 



FIGURE 8- PLATFORM RESPONSE DURING 
ANGLE REGULATION 

Angle Regulation Simulation: DC Bus Response 



time (sec) 


FIGURE 9- DC BUS DURING ANGLE 
REGULATION 
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system is in discharge mode, with a minimal electrical 
load. Simulation parameters are summarized in row 2 
of Table 2. 

The attitude regulator (see Fig. 3) is a simple 
PID, with gains set heuristically to limit the flywheel 
motor currents to 30 A. In this simulation, a one 
degree step change is made in theta *, which is the 
angle commanded into the attitude regulator. The 
results are summarized in Figure 8. Note that the 
angle settles to the set point and the torque settles 
back to zero within a few seconds; this represents 
satisfactory response to a step change in commanded 
angle. 

The bus regulation during this step change in 
angle is displayed in Figure 9. Since the system is in 
discharge mode, all of the current for the load is 
provided by the flywheel system. In deriving equation 
(11), it is assumed that the q-axis rotor reference 
frame current, f ?s , is in steady state; however, with a 
step change in angle, the q-axis current jumps 
abruptly, causing some error in the DC current. There 
is a brief spike in the I DC current, and a resultant dip in 
bus voltage. Note that the system responds very 
quickly to the error; the bus voltage is regulated to 
within 0.5 volt of the commanded value, which is less 
than 0.5% of the nominal bus voltage setpoint. 


x iq* Angle Regulation Simulation: Flywheel Rotor Speeds 



time (sec) 


FIGURE 10- FLYWHEEL SPEEDS DURING 
ANGLE REGULATION 



FIGURE 11- PLATFORM SPEED DURING 
ANGLE REGULATION 
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Torque Disturbance Rejection 


A physical understanding of the dynamics of the 
platform motion can be achieved by examining the 
flywheel speeds and the platform torque and angular 
velocity during the step change. A plot of the speeds 
of the two flywheels during the step in commanded 
angle is shown in Figure 10. In this plot, the HSS 
flywheel speed is plotted as the negative of its actual 
speed, so the two speeds can be observed on the 
same plot. The same torque command results in a 
greater speed change in the HSS as compared to D1 
due to the differences in the flywheels’ inertia. From 
about 0.5 to 1.3 seconds, the HSS flywheel 
accelerates in the negative speed direction and the 
D1 flywheel decelerates, generating a positive 
platform torque and a positive angular velocity of the 
platform. After 1.3 seconds, the flywheel 
accelerations change direction, exerting a net 
negative torque on the platform, slowing it down as it 
reaches the position setpoint. 

Figure 11 shows the corresponding platform 
angular velocity and the total torque, T p . 

Torque Disturbance Rejection 

The final simulation demonstrates the 

disturbance rejection capability of the controller. A 
step change in torque is applied to the platform 
dynamics (the "T d " input to the "1 DOF Platform Plant 
Model" of Fig. 3) and the capability of the controller to 
maintain the platform angle at the setpoint is studied. 
The system is in discharge mode, with a minimal 
electrical load. Simulation parameters are 

summarized in row 3 of Table 2. This represents a 
physical disturbance to the spacecraft. 

The applied torque disturbance is 0.05 N*m 
which is about 10% of the maximum torque capability 
of either flywheel. Figure 12 shows the application of 
the torque disturbance at 0.5 seconds. The sum of 
the torque on the platform thus jumps from zero to 
0.05 N*m at 0.5 seconds. This causes the platform 
to begin to rotate away from its setpoint of 0 degrees. 
The resulting position error is acted upon by the PID 
in the attitude regulator block (figure 3) and creates 
the appropriate current commands, f r qsI and I* r qs2 , to 
the flywheels as shown in Fig. 13. This results in a 
net torque applied to the platform to counteract the 
disturbance torque and reduce the position error as 
shown in Fig. 12. 

Figure 14 shows the speed of each of the 
flywheels as the disturbance torque is rejected and 
the platform position is returned to its original position. 

CONCLUSIONS 

This paper presents the simulation results of a 
combined flywheel energy storage system and a 
single axis attitude control system. The simulation 
demonstrates that it is possible to use two flywheels 
to simultaneously provide energy storage for the 
electrical system while regulating the spacecraft 
platform attitude in a single axis. A new control 
algorithm was developed which uses the platform 
angle requirement and the electrical system power 
requirement to determine the appropriate current 



FIGURE 12- PLATFORM DURING TORQUE 
DISTURBANCE REJECTION 



Time (Seconds) 

FIGURE 13- DC BUS DURING TORQUE 
DISTURBANCE REJECTION 



FIGURE 14- FLYWHEEL SPEEDS DURING 
TORQUE DISTURBANCE REJECTION 


commands to the two flywheel motor controllers. A 
complete system model was developed to test the 
new controller in simulation. The model includes both 
the electrical and mechanical components: the 
flywheel motor/generator units and their control, the 
DC bus and the electrical load, the platform inertial 
model and the new control algorithm. Simulation 
results show that the new controller can successfully 
maintain commanded platform angle through all mode 


NASA/TM— 2002-21 1812 


7 


transitions, and can regulate the DC bus voltage 
under charge reduction and discharge conditions 
while maintaining the correct platform angle. In 
addition, a step change in platform angle can be 
achieved with minimal impact (less than 0.5%) on the 
DC bus voltage. The system also can reject step 
torque disturbances. 

Future work includes demonstrating this control 
algorithm experimentally using two flywheel units, an 
air table and an electrical load. Additionally, more 
work will be done to optimize the controller to achieve 
the best possible overall system response. 

This work is a first step toward an overall 
program which will ultimately reduce the cost and 
mass of future space vehicles with the combination of 
the energy storage and the attitude control systems. 
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